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Abstract 



A nonrigid rotor model is developed from the two-parameter quantum algebra U qp (u2)- [This 
model presents the U qp (u2) symmetry and shall be referred to as the gp-rotor model.] A ro- 
tational energy formula as well as a gp-deformation of E2 reduced transition probabilities are 
derived. The gp-rotor model is applied (through fitting procedures) to twenty rotational bands 
of superdeformed nuclei in the A ~ 130, 150 and 190 mass regions. Systematic comparisons 
between the gp-rotor model and the g-rotor model of Raychev, Roussev and Smirnov, on one 
hand, and a basic three-parameter model, on the other hand, are performed on energy spec- 
tra, on dynamical moments of inertia and on B(E2) values. The physical signification of the 
deformation parameters q and p is discussed. 
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1 Introduction 



Quantum groups and quantum algebras, introduced at the beginning of the eightees, 1-5 continue 
to attract much attention both in mathematics and physics. For the Physicist, a quantum 
algebra is commonly considered as a deformation (g-deformation) of a given Lie algebra. During 
the last four years, several works have been performed on two-parameter quantum algebras and 
quantum groups (gp-deformations). 6 ~ 17 (For an elementary introduction to q- and gp-quantum 
algebras, the reader should consult Ref. 14.) 

Most of the physical applications, ranging from chemical physics to particle physics, have 
been mainly concerned up to now with one-parameter quantum algebras (g-deformations). 
In particular, in nuclear physics we may mention applications to rotational spectroscopy of 
deformed and superdeformed nuclei, 18-26 to the interacting boson model, 27,28 to the Moszkowski 
model, 29 ' 30 to the U(3) shell model 31 and to the Lipkin-Meshkov-Glick model. 32 There exist also 
applications to particle physics, as for example to quote a few, to hadron mass formulas 33 ' 34 and 
to Veneziano amplitudes. 35 ' 36 Among the just mentioned applications, only the ones in Refs. 25 
and 36 rely on the use of two-parameter deformations. 

The aim of the present paper is two-fold: (i) to further develop and (ii) to apply (to rotational 
bands of superdeformed nuclei) the nonrigid rotor model briefly introduced in Ref. 25. The 
latter model, referred to as the gp-rotor model, is based on the two-parameter quantum algebra 
U qp (u 2 ) while the g-rotor models introduced by Iwao 18 and Raychev, Roussev and Smirnov 19 (see 
also Refs. 20-23 and 26) are based on the one-parameter quantum algebra U q (su 2 ). One of the 
objectives of this work is to show what we gain when introducing a second "quantum algebra" - 
type parameter, i.e., when passing from the U q (su2) symmetry to the U qp (u2) symmetry. 

The organization of this paper is as follows. The gp-rotor model is introduced in Sec. 2. 
Subsection 2.1 deals with the mathematical ingredients of the model. The gp-rotor model 
itself is developed in Subsec. 2.2 (rotational energy formula) and in Subsec. 2.3 (E2 transition 
probabilities). Section 3 is devoted to the application of the gp-rotor model to the description 
of superdeformed (SD) bands of nuclei in the A ~ 130, 150 and 190 mass regions. The results 
obtained from the gp-rotor model for rotational energy spectra, dynamical moments of inertia 
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and B(E2) values are compared to the ones derived from the g-rotor model and from a basic 
(d la Bohr-Mottelson) model. Finally, some concluding remarks are presented in Sec. 4. 



2 A gp-Rotor Model 

2.1 The quantum algebra U qp (u2) 

The quantum algebra U qp (u2) can be constructed from two pairs, say {d+,d + } and {dt,d_}, 
of gp-deformed (creation and annihilation) boson operators. The action of these gp-bosons on 
a nondeformed two-particle Fock space {\n + ,n_) : n + G N, n_ G N} is controlled by 



a+ \n+, n_) = \j [[ra+ + - + h] qp \n+ + l,n. 



a+\n+,n_) = J[[n+ + - - ^}} gp \n+ - 1, ra_), 



dl|n+,n_) = J[[n- + ^ + ^]] qp |n+,n_ + 1), 



d_ |n+, n_) = J [[n_ + - - ^}} qp \n+, n_ - 1). (1) 



In the present paper, we use the notations 



and 



[MW : = "-^zy ( 2 ) 

== [[*]]«- = ( 3 ) 
where X may stand for an operator or a (real) number. For Hermitean conjugation require- 
ments, the values of the parameters q and p must be restricted to some domains that can be 
classified as follows: (i) q G R and p G R, (ii) q G C and p G C with p = q* (the * indicates 
complex conjugation), and (iii) q = p -1 = e 1/3 with < (3 < 2n. The two pairs {d+,d + } and 
{dl,d_} of gp-bosons commute and satisfy 

a±a+ = [[N± + l]] qp , d+d± = [[N±\] qp , (4) 



where N + and N_ are the usual number operators with 

N± |n + ,n_) = n± |n+,n_). (5) 

Of course, the gp-bosons a± and a± reduce to ordinary bosons (denoted as alj: and a± in Refs. 37 
and 38 and in Subsec. 2.3) in the limiting situation where p = q^ 1 — > 1. 

The passage from the (harmonic oscillator) state vectors |n+,n_) to angular momentum 
state vectors \I,M) is achieved through the relations 

I:=^(n+ + n_), M:=^(n + -n_) (6) 

and 

\I,M) = \I + M,I-M) = |n+,n_). (7) 
Equations (|T|) may thus be rewritten as 



4 \I,M) = J[[I±M+± + ±]] qp \l+l,M±± 



s± |J,M) = y[[/±M + i-i]] w |/-~ M T ^), (8) 

so that the gp-bosons behave as ladder operators for the quantum numbers / and M (with 
\M\<I). 

We are now in a position to introduce a gp-deformation of the Lie algebra u 2 . A simple 
calculation shows that the four operators J a (a = 0, 3, +, — ) given by 

Jo := ^(N+ + N.), J 3 := ±(N + -N-), J + : = J_ := ~a + _~a + (9) 

satisfy the following commutation relations 14 ' 39 

[J 3 ,J±] = ±J±, [J+,J_] = (gp) J °- j3 [[2J 3 ]W, [Jo,Ja] = 0. (10) 

We refer to t/ gp (u 2 ) the (quantum) algebra described by fllP]) . To endow t/ gp (u 2 ) with a Hopf 
algebraic structure, it is necessary to introduce a co-product A qp . The latter co-product is such 



that: 



\p( Jo) = Jo ® 1 + 1 ® Jo, 
A W (J 3 ) = J 3 ® 1 + 1 ® J 3 , 

A gp (J ± ) = J ± ® (gp)^ Jo (gp- 1 ) +|Ja + (gp)^ Jo (sp" 1 )"*' 8 ® J± (11) 



and is clearly seen to depend on the two parameters q and p. [Note that with the constraint 
p = q* , to be used in Subsec. 2.2, the co-product satisfies the Hermitean conjugation property 
(Agp(J±))T = A pq (J T ) and is compatible with the commutation relations for the four opera- 
tors A qp (J a ) (with a = 0,3,+,—).] The universal 7£-matrix (for the coupling of two angular 
momenta I = \) associated to the co-product A qp reads 



1Z 



(p 0^ 

y/pq 

p — q y/pq 

\0 pj 



(12) 



and it can be proved that TZ pq verifies the so-called Yang-Baxter equation. 
The operator defined by 14 

C 2 (U qp (u 2 )) := l -(J + J_ + J_J + ) + I [[2}] qp (qp) J °- J * ({[J 3 ]] qp ) 2 (13) 

is an invariant of the quantum algebra f/ gp (u 2 ). It depends truly on the two parameters q 
and p. The invariant C2{U qp (u2)) will be one of the main mathematical ingredients for the qp- 
rotor model to be developed below. Hence, it is worth to examine its structure more precisely, 
especially its dependence on two independent parameters. Equation QTTD suggests the following 
change of parameters 

Q:=(<7P~¥, P:=(gp)i (14) 
Then, by introducing the generators A a (a = 0, 3, +, — ) 

A := J , A 3 := J 3 , A ± := {qp)-^ J ^ J ± , (15) 

it can be shown that the two-parameter quantum algebra ?7 ? p(u 2 ) is isomorphic to the central 
extension 

U qp (u 2 ) = ui <g> Uq(bu 2 ), (16) 

where Ui is spanned by the operator A and Uq(su 2 ) by the set {A 3 , A + , A_}. The Q- 
deformation Uq(su 2 ) (a one-parameter deformation!) of the Lie algebra su 2 corresponds to 
the usual commutation relations 

[A 3 ,A ± ] = ±A±, [A + ,A_] = [2A 3 ] Q . (17) 
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Furthermore, the co-product relations ([LTD leads to 

A 9P (J ± ) = P A ^o)-h Ag(i4±)) (1 8) 

where the co-product Ag is given via 

A Q (A>) = A ® 1 + 1 ® A), 
A Q (A) = A ® 1 + 1 ® A, 

A (At) = At ® Q +Aa +g^ 3 <g> At. (19) 

Equations (|T7| ) involve only one parameter, i.e., the parameter Q. However, two parameters 
(Q and P) occur in (|T8| ) as well as in the invariant (^(L/^u^) transcribed in terms of Q and 
P. As a matter of fact, (TTSf) can be rewritten as 



C 2 (U qp (u 2 )) = P 2A ^ 1 C 2 (C/ Q (su 2 )), (20) 

where 

C 2 (C/ Q (su 2 )) := ~ (A, A_ + A_A + ) + ~ [2] Q ([A]q) 2 (21) 

is an invariant of ?7q(su 2 ) [compare Eqs. ([13]) and (|2T|)] . As a consequence, of central importance 
for the gp-rotor model of Subsec. 2.2, the invariant C 2 (t/ gp (u 2 )), in either the form ( [T3|) or the 



form (pD]) , depends on two parameters. Finally, it should be noted that C 2 ([/q P (u 2 )) can be 
identified to the invariant of L/g(su 2 ) and to the Casimir of su 2 when p = q~ l and p = q~ l — > 1, 
respectively. In this sense, the ^7 qp (u 2 ) symmetry encompasses the C/ g (su 2 ) and su 2 symmetries. 

To close this section, let us say a few words on the representation theory of £/q P (u 2 ) in the 
case where neither q nor p are roots of unity. An irreducible representation of this quantum 
algebra is described by a Young pattern [<pi\ <£> 2 ] with (pi — <^ 2 = 2/, where 2/ is a nonnegative 
integer (J will represent a spin angular momentum in what follows). We note \[(pi', <£> 2 ], M) 
(with M = —I, —I + 1, • • ■ , +/) the basis vectors for the representation [(pi, (p 2 ]. 

We are now ready to develop a gp-rotor model for describing energy levels and transition 
probabilities for deformed and superdeformed nuclei. 
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2.2 Energy levels 

We want to construct a nonrigid rotor model. As a first basic hypothesis (Hypothesis 1), we 
take a rigid rotor with U qp (u2) symmetry, thus introducing the nonrigidity through the qp- 
deformation of the Lie algebra u 2 - More precisely, we assume that the qp-rotor Hamiltonian H 
is a linear function of the invariant C 2 (U qp (u 2 )): 

H = ^C 2 {U qp (u 2 )) + E , (22) 

where X denotes the moment of inertia of the rotor and Eq the bandhead energy. As a second 
hypothesis (Hypothesis 2), we take (pi = 21 and (p 2 = 0. This means that we work with state 
vectors of the type \I,M) = \[2I;0], M). Therefore, the eigenvalues of H are obtained by 
the action of H on the physical subspace {|J, M) : M = —I, —I + 1, • • • , +/} of constant 
angular momentum I corresponding to the irreducible representation [21; 0] of U qp (u 2 ). The 
two preceding hypotheses lead to the energy formula 

E(I) qp = ±[[I]] qp [[I + l]] qp + E (23) 

for the gp-deformed rotational level of angular momentum /. 
By introducing s = \nq and r = hap, Eq. (^) yields 

v )ip 2 j sinh 2 (^) V ; 

Preliminary studies have lead us to the conclusion that a good agreement between theory and 
experiment cannot be always obtained by varying the parameters s and r (or q and p) on the 
real line R, a fact that confirms a similar conclusion reached in Ref. 20 for p = q^ 1 G R. In 
addition, if we want that our gp-rotor model reduces to the g-rotor model developed by Raychev, 
Roussev and Smirnov 19 when p = q^ 1 (or equivalently r = —s), we are naturally left to impose 
that (s + r) and (s — r)/i should be real numbers. [Observe that the two constraints (s + r) G R 
and (s + r)/i G R ensure that the energy E(I) qp is real as it should be.] Furthermore, we 
shall see that for certain SD bands, a good agreement between theory and experiment requires 
that the parameters s and r vary on the real line R. Thus, we shall consider the two possible 
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parametrizations: 



S ~\~ T S — V 

— - — = P cos 7 G R, = (3 sin 7 G R, 

_ _ 1 



(25) 



(b) = /3cos7GR, = G 1R, 

\ ) 2 r ' 2i i 



so that the parameters g and p read 



g — e /3cos7 e +i/3sin7 p = q* — g/ 3 ™ 8 "? g-i/Ssin7 



g — g/?cos7 e +/3sin7 ^ _ g/3 cos 7 e -/?sin7 



(26) 



Thus, our gp-rotor model involves two independent real parameters (3 and 7 corresponding 
either to (a) the two complex parameters q and p subjected to the constraint p = q* or to (b) 
the two real parameters q and p. In terms of the parameters (3 and 7, the energy formula (p4j) 
takes the form 



pm 1 J2/-D/3COS7 sin(J/3sin 7 ) sin[(I + l)/3sin 7 ] 

21 sin (p sin 7) 



or 



E{1) = J_ e ( 2 /-i )/3 cos7 sinh(^sin 7 ) sinh[(J + l)/?sin 7 ] + ^ 
2X sinh (/? sin 7) 

in the parametrizations of type (a) or (b), respectively. We shall use both Eqs. (27a) and (27b) 
in our fitting procedures. 

In the (a)-parametrization, to better understand the connection between our gp-rotor model 
and the g-rotor model of Ref. 19, we can perform a series analysis of Eq. (27a). A straightforward 
calculation allows us to rewrite Eq. (27a) as 

1/00 CO \ 

E{I) m * = TyT E 7) [/(/ + l)] n + (2/ + 1) E Cntf, 7) [/(/ + 1)]" + Eo, (28) 

Z1 /37 \n=0 n=0 / 

where the expansion coefficients c? n (/3, 7) and c n (P,j) are given in turn by the series 

o2n-l 



sin {(3 sin 7) 



(2k + 2n)\ k\n\ 



>2n-l 



x £ {(cos 7 ) 2fc+2n cos(/5sin 7 ) - cos[(2£; + 2n) 7 ]|- ^ + ^ 

t {(cos,)— cos(/*sin 7 ) - cos[(2* + 2„ + l) T ]} , 29) 



X 

A;=0 



sin 2 (/3 sin 7) 

^ r 1 ^2fc+2n+l (jfe + n )! 

(2fc + 2n+l)! k\n\ 



In Eq. (|28|) , we have introduced the deformed moment of inertia: 

X M =Te 2/3cos7 , (30) 

which gives back the ordinary moment of inertia when 7 = tt/2 (i.e., q = p^ 1 = e 1/3 ). In the 
limiting situation where 7 = tt/2, the coefficients c n (/3, 7) vanish and the energy formula ([28|) 
simplifies to 

1 00 on— 1 

-^-i- Sill ^J? „ ^ Tt. 



where j n -\ denotes a spherical Bessel function of the first kind. Equation ( pl|) was derived 
by Bonatsos et at 20 for the g-rotor model with q = e 1/3 in order to prove the mathematical 
parentage between the g-rotor model and the variable moment of inertia (VMI) model. 40-42 
The series (|3lD corresponds indeed to the compact expression 

E(I) qq -. = ± [I] q [I + l] q + E , (32) 

to be compared with Eq. fl23|). Note that Eq. (32) corresponds also to the (b)-parametrization 
with 7 = 7r/2. The transition from Eq. (|23"D to Eq. ( p2[) illustrates the descent from the U qp (\i2) 



symmetry of the gp-rotor to the U q (su>i) symmetry of the g-rotor. A further descent in symmetry 
is obtained when (3 — > (i.e., q = p~ x — > 1): in this case [I) q [I + l) q — > 1(1 + 1) and we get 
[from Eq. (j32|) 1 the usual energy formula corresponding to the rigid rotor with SU2 symmetry. 

2.3 E2 transition probabilities 

We now examine the implication of the U qp (xi2) symmetry on the calculation of the electric 
quadrupole transition probability T(E2; / + 2 — > I). Let us start with the ordinary expression 



9 



of the reduced transition probability, namely, 



B(E2;I + 2 -> I) 



16tt 



Ql 



7 + 2,M,2,0|/ + 2,2, J,M) 



(33) 



for an E2 transition. 43 In Eq. (|33|) , Qo is the intrinsic electric quadrupole moment in the body- 
fixed frame. The coefficient of type (j,m,k, fi\j,k,j',m') in the right-hand side of Eq. (j33|) 
is a usual Clebsch-Gordan coefficient for the group SU(2). Our goal is to find a gp-analog of 
T(E2; 1+2 — > I) and, thus, of Eq. (|33"|). The strategy for obtaining a qp- analog of _B(E2; 1+2 — > 
I) is the following: 

(i) We first rewrite the SU(2) Clebsch-Gordan coefficient of Eq. fl33|) in terms of a matrix 
element of an SU(2) unit tensor operator t kfMa with k — 2, p, = and a = —2. This may be 
done from the general formula 38 



(f, m'\t ktMX \j, m) = 5(f, j + a)5(m', m + /i)(-i) 2fc (2j" + (j, m, k, n\j, k, f, m'), (34) 

which shows that the irreducible tensor operator tk^a produces the (angular momentum) state 
vector \j + a,m + fi) when acting upon the state vector \j, m). Then, Eq. (|33|) is amenable to 
the form 



B(E2;I + 2 -> I) 
by making use of Eq. 



16n 



Ql (2/+1) 



(I,M\t 20 „ 2 \I + 2,M) 



(35) 



(ii) We know that the general operator tk^a can be realized in terms of two pairs {a^, a + } 
and {al, a_} of ordinary boson operators. In this respect, we may consider the so-called van der 
Waerden 38 realization of tk m - There are several ways to gp-deform the operator tk m - Here, we 
choose to define a qp- deformation tk^qp) by replacing, in the van der Waerden realization of 
tk/j.a, the ordinary bosons {ajl, a±} by gp-deformed bosons {ajl, a±} and the ordinary factorials 
x\ by gp-deformed factorials = [[x]] qp [[x — 1]] 



qp 



qp 



for x G N. We thus obtain 



t kf u*{qp) = (-1)* 



-a{[[k + »}W- [[k-fAW [[k + a}} qp \ [[k-a]] qp \ [[2j - k + a]] qp \\ I 



E(-i) 



[[2j + k + a + \}} qp \ 
(a±) k -^- z {a_) k - a - z {aXY +a+z (h + ) z 
[[k-fi- z]]gp\ [[k — a — z]]gp\ [[n + a + z]] qp \ [[z}\ 



(36) 



qp ■ 



10 



In particular, the gp-deformed operator £20-2(9^) connecting the state vector \I + 2,M), with 
j = / + 2, to the state vector \I, M), with f = I, reads 



^^Wm^) i(W - (37) 

an expression of direct interest for deriving the gp-analog of B(E2; I + 2 — > I). 

(hi) We assume that the qp- analog B(E2; I + 2 — > I) qp of S(E2; / + 2 — > /) is simply 
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B(E2;I + 2^I) qp := — Q 2 [[21 + 1]} 



qp 



(I,M\t W - 2 (qp)\I + 2,M) 



2 



(3? 



167T 

[Equation (^) constitutes the third and last hypothesis (Hypothesis 3) for our gp-rotor model.] 
By using Eqs. ([37]) and (§), the relevant matrix element of the operator £20-2 (<?£>) is easily found 
to be 

(I,M\t 20 ^(qp)\I + 2,M) = 

• [[3}} qp [[A]] qp [[2I]] qp \ [[I + M+ l]] qp [[I-M+ l]] qp [[I + M + 2}} qp [[I - M + 2}} qp 



[[2/ + 5]] gp \ 



(39) 

Then, the introduction of Eq. (|39"D into Eq. (|38D yields 



R(F2T , 2 n 5 2 [[3}} qp [[4}) qp ([[/+i]W[[/ + 2]W) 2 

1 ' >qp 16^° [[2]] qp [[2I + 2]} qp [[2I + 3}} qp [[2I + 4]] qp [[2I + 5}} qp { ] 

in the case of the K = M = bands. 

For the purpose of comparison with experimental results, we must calculate the E2 transition 
probability T(E2; I + 2 I) in the gp-deformed scheme. We define such a probability by 

T(E2; I + 2 -> I) qp := 1.223 10 9 (£ 7 (I + 2) w ) 5 5(E2; J + 2 -> (41) 

Equation ( |4T| ) turns out to be a simple gp-deformation of the usual E2 transition probability. 
[In Eq. ©, T(E2;/ + 2 ->• is in units of sec" 1 , S 7 (/ + 2) ?p := E(I + 2) qp - E(I) qp in units 
of MeV, and 5(E2; / + 2 -> in units of e 2 fm 4 .] 

At this stage, a contact with the formula -B(E2; 1+2 — > I) q derived by Raychev, Roussev and 
Smirnov 19 is in order. First, by taking p = g _1 the right-hand side of (|40|) may be specialized 
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to the expression of B(E2; I + 2 — > I) q obtained in Ref. 19. Hence, our gp-rotor model for the 
E2 transition probability admits as a particular case the corresponding g-rotor model worked 
out in Ref. 19. Second, it can be shown that 

£?(E2; I + 2 - I) qp = P~ A ^B(E2; I + 2 - I) Q , (42) 

where P and Q are given by (PI). 

Let us close with a remark. Should we have chosen to find a gp-analog of the Clebsch-Gordan 
coefficient in (|33|), we would have obtained 39 



(/ + 2, M, 2, 0\I + 2, 2, J, M) qp = (1 + 2, M, 2, 0| J + 2, 2, /, M) Q (43) 
and, consequently 

B(E2; I + 2^I) qp = B(E2; 1 + 2^ I) Q . (44) 
We prefer to use (|42|) rather than ( 44]) because the factorization in (|42| ) parallels the one in 



3 Description of Superdeformed Bands 
3.1 Fitting procedure 

The gp-rotor model developed in Sec. 2 was applied to twenty rotational SD bands of nuclei in 
the A ~ 130, 150 and 190 mass regions. The 7-ray energies 

£ 7 (J) := E(I) - E(I - 2) (45) 

were computed from the energy formula 

E(I) . E(I) qq . = i e<-»« M ""^ + m + So (46a) 



or 



£(/) = £(/) ?p - - e + £ 



(46b) 



that correspond to Eq. (27a) or (27b), respectively, with 

a := /5cos7, b := /3sin7. (47) 
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The free parameters of the gp-rotor model are then a, b and X. 

For the sake of comparison, we also computed the transition energies E y (I) and performed 
an analysis of the same SD bands from two other models. First, we used 

E(I) = E{I) qq -, = — + E , (48) 

where f3' (defined by q = e 1/3 ') and X' are the two free parameters of the g-rotor model of Ref. 19. 
Second, we also applied the energy formula 

E(I) = A' 1(1 + 1) + £'[/(/ +1)] 2 + C'[I(I+1)] 3 + E (49) 

arising from the (Bohr-Mottelson) basic model 43 restricted to three free parameters, i.e., A', B' 
and C 

For each of the three models, the parameters (a, b and X for the gp-rotor; j3' and X' for the 
g-rotor; A', B' and C for the basic model) were fixed by minimizing 



X 



— E 

\ n — m j 



E?(I)-E?(I) 
A£ 7 (J) 



2 



(50) 



where n is the number of experimental points included in the fitting procedure, m is the number 
of freely varied parameters, and AE y (I) are the experimental errors. 

3.2 Results and discussions 
3.2.1 Fitting of data 

We present in Table 1 the free parameters and the x-values, for the q- and the gp-rotor models, 
obtained from the twenty fitted SD bands. (The tables and figures of this paper can be obtained 
from the authors.) For space saving purposes, we do not report the corresponding results 
obtained with the basic model since the x- va l ues are generally higher than the ones derived 
from the gp-rotor model (a fact to be confirmed in Subsec. 3.2.2). 

Table 1 exhibits two general trends: first, the best results are obtained in the A ~ 190 mass 
region for the two models and, second, the x-values for the gp-rotor are better than those for 
the g-rotor. Indeed, the x _va l ues obtained for the gp-rotor (respectively, g-rotor) are between 
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0.6 and 4.2 (respectively, 0.9 and 29.4) in the A ~ 190 mass region except for 194 Hg(a) with 
X = 9.1 (respectively, 35.2) while the x _va hies are between 1.9 and 20.3 (respectively, 8.8 and 
87.9) in the two other mass regions. The high values obtained for \ are not surprising: in 
the standard definition of \i Eq. (|50|), the difference between the theoretical and experimental 
transition energies is divided by the experimental error AE y (I) that is equal to 0.5 keV except 
for the recent experimental data 44,45 on 192 Hg, 194 Hg(a) and 194 Hg(b) for which Ai? 7 (J) is as 
below as 0.1 keV. Therefore, we may emphasize the excellent quality of the fits for the A ~ 190 
bands, especially in the case of the gp-rotor model. 

For the gp-rotor model, the quality of the fits is connected to the nature (real or complex) 
of the parameters q and p. The best fits were obtained by taking: (i) the (a)-parametrization 
(i.e., q = e a+lb and p = e a ~ lb ) for 146 Gd and the 190 SD bands and (ii) the (b)-parametrization 
(i.e., q = e a+b and p = e a ~ b ) for the 130 and 150 SD bands. To illustrate our results, we globally 
characterize in Fig. 1 the twenty SD bands by their position in the plane of the two "quantum 
algebra" -type real parameters a and b. As it was shown in Ref. 20, the parameter f3' of the q- 
rotor (which occurs in a sine like the parameter b in the (a)-parametrization of the gp-rotor) can 
be interpreted as a softness or stretching parameter of the nucleus, similar to the parameter 
a of the VMI model. 40-42 We adopt this interpretation for the parameter b (that coincides 
with the parameter (3' when a = 0) in the (a)-parametrization. Then, the (b)-parametrization 
describes a distortion phenomenon (decrease of the dynamical moment of inertia with the spin 
of the nucleus) rather than a stretching phenomenon. In the (a)- and (b)-parametrizations, the 



role played by the parameter a, appearing in the exponential term of Eq. ([46]), is clear. The 
parameter a has a crucial effect of correction on the distortion (stretching or anti-stretching) 
function of the parameter b. In addition, we note from Table 1 that the sign of a is the same as 
that of the difference X — X' . In other words, at high angular momenta, the exponential term in 
(|46|) moderates (when a < 0) or accentuates (when a > 0) the contribution to the energy of ^ 
with respect to the contribution of ^7. Therefore, the parameter a can moderate or accentuate 
the distortion phenomenon of the nucleus. 

Before performing a systematic comparison between the three models under consideration, 
we present in Tables 2-6 the calculated and experimental transition energies for the gp-rotor 
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model. The numerical results in Tables 2-6 confirm the preceding interpretation of the free 
parameters of the gp-rotor model. Here again, we note that the quality of the fits is better 
in the A ~ 190 region than in the A ~ 130 and 150 regions. This reflects the fact that the 
7-ray energies range from 200 to 900 keV (respectively, 600 to 1700 keV) for angular momenta 
ranging from 8 to 50 (respectively, 14 to 64) for A ~ 190 (respectively, A ~ 130 and 150). 

3.2.2 Comparative analyses 

In order to confirm the difference (already evocated in the x _vames analysis) between the 
gp-rotor model and the basic model, we consider three representative nuclei for each of the con- 
sidered mass regions. Figure 2 shows the differences between the calculated and experimental 
transition energies for the nuclei 132 Ce, 152 Dy and 192 Hg obtained from the basic and gp-rotor 
models. It is clear that the gp-rotor model is more appropriate, in particular for 132 Ce, for 
describing the distortion phenomenon than the basic model. Therefore, we switch to a detailed 
comparison between the q- and gp-rotor models. Figures 3-7 display the results (in terms of 
differences as in Fig. 2) afforded by the q- and gp-rotor models for the twenty SD bands un- 
der study. Two remarks arise from Figs. 3-7. First, the preceding x-values analysis is clearly 
confirmed. Second, we observe that the gp-rotor model is much better than the g-rotor one 
when the distortion phenomenon is particularly pronounced. For example, in the case of the 
192 Hg band that presents nineteen transitions and where the variation of the stretching effect 
becomes less important at high spin, the gp-rotor model provides the best results. 

An alternative way to analyse the stretching phenomenon in the A ~ 190 region amounts 
to compare the theoretical and experimental dynamical moments of inertia 



theoretical energies E in ([51]) as given by (|46| ) (respectively, (f48|)) for the qp- rotor (respectively, 




(51) 



and 




g-rotor) model. [The dynamical moments in floTf) and [5^) are in units of h 2 MeV .] Figure 8 
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shows the results for four SD bands of the three nuclei 19 °- 192 - 194 Hg: the experimental moments 
of inertia are calculated from Refs. 44, 45 and 54 and the theoretical ones by using the free 
parameters ~f, a, b and ^7, f3' for the U qp (u2) and U q {su2) symmetries, respectively. The 
gp-rotor results are much closer to the experimental results than the g-rotor ones, due to the 
influence of the parameter a. [In passing, Fig. 8 shows that globally, both for the q- and qp- 
rotor models, the second derivative of the energy is significative when calculated with the fitted 
values of the free parameters.] 

A last way to compare the gp-rotor model with the two others is to use experimental 
values of E2 transition probabilities. From such values, we can compute two different intrinsic 
electric quadrupole moments, namely, (Qo) qp and (Qo) q for the t/ gp (u 2 ) and ?7 g (su 2 ) symmetries, 
respectively. For the U qp (u2) symmetry, (Qo) qp is deduced from (|40|) and (|4]), where we take 
the experimental value for the E2 transition probability and all the other terms (including the 
transition energies) are calculated from the parameters of the gp-rotor model obtained from the 
optimization of energy. A similar calculation is conducted for (Qo) q corresponding to the ^7 9 (su 2 ) 
symmetry. The experimental intrinsic electric quadrupole moment (Qo)cx corresponds to the 
su 2 symmetry: it is calculated from ( (4"0"| ) and ( f4"T|) with q = p~ l — > 1 by taking the experimental 
E2 transition probability and the experimental 7-ray energy. We present in Tables 7-9 the 
values of the quadrupole moments (Qo) qp , [Qo)q an d (Qo)ex, together with the experimental 
errors, computed for 192 Hg and 194_196 Pb with the experimental E2 transition probabilities 
of Refs. 57, 60 and 61. We see that the values of {Qo) qp are in better agreement with the 
experimental quadrupole moments (Qq) cx than the values of {Qo) q - To further compare the 
symmetries £/q P (u 2 ), ?7 g (su 2 ) and su 2 , it is interesting to calculate the geometrical factor of type 

G[I) ;= l ^(E2;/+2^) m 

for the three symmetries. Figure 9 displays this factor as a function of the spin / for the 
three nuclei 192 Hg and 194 - 196 pb. At high spin, the increasing of G(I) characterises the two 
"quantum algebra" -type models, while G(I) reaches a limit value for the su 2 symmetry. Note 
that G(I) increases less strongly (i.e., more linearly) for the t/ gp (u 2 ) symmetry than for the 
[/q(su 2 ) symmetry when the parameter b differs from the parameter (3' . 
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4 Conclusions 



In this paper, we concentrated on a new nonrigid rotor model (the gp-rotor model) based on 
three hypotheses in the framework of an investigation of the two-parameter quantum algebra 
Uqp(\i2). The two facets of this model consist of a three-parameter energy level formula and 
a gp-deformed E2 transition probability formula. As limiting cases, the gp-rotor model gives 
back the g-rotor model 19 (when p = g _1 ) based on the quantum algebra U q (su 2 ) and the rigid 
rotor model (when p = g -1 — > 1) based on the Lie algebra su 2 . 

Twenty rotational bands of superdeformed nuclei in the A ~ 130, 150 and 190 mass regions 
were used to test our gp-rotor model and to compare it to the g-rotor model and to a basic 
(with a three-term polynomial energy formula) model. The main results may be summarized 
as follows. First, the gp-rotor model is better than the g-rotor model and the basic model as far 
as energy spectra are concerned. Second, the energy fits for the twenty SD bands are in good 
agreement with experiment both for the q- and gp-rotor models. However, a marked difference 
between the latter two models manifests itself in the energy spectrum and also in the second 
derivative of the energy (i.e., for the dynamical moment of inertia). Third, in terms of B(E2) 
values the results afforded by the U qp (u 2 ) symmetry are between those given by the C/q(su 2 ) 
symmetry and the su 2 symmetry: the B(E2) values for the qp- rotor model increase more or 
less linearly with spin, a result that does not hold for the g-rotor model. 

As a general conclusion, the gp-rotor is appropriate for describing the collective phenomenon 
of distortion occurring in the rotation of the nucleus (increase or decrease of the dynamical 
moment of inertia with the spin). The net difference between the q- and gp-rotor models 
comes from the "quantum algebra" -type parameter a that tends to smooth the (spherical 
or hyperbolical) sine term in the energy and thus accentuates or moderates the distortion 
phenomenon of the nucleus. 

To close this paper, let us mention that Hypothesis 2 (i.e., (p 1 = 21 and (p 2 = 0) of our 
model might be abandoned. This would lead to a a la Dunham formulation for describing more 
complicated rotational spectra of deformed and superdeformed nuclei or rovibrational spectra 
of diatomic molecules. As a further extension, it would be also interesting to combine our model 
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with one of Ref. 24 (based on the g-Poincare symmetry) in the case of heavy nuclei. Work in 
these directions is in progress. 
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Table captions: 

Table 1. Free parameters for the q- and gp-rotor models: f3' corresponds to q — e 1/3 '; a = (3 cos 7 
and 6 = /9sin7 correspond: (i) to q — e a+lb and p = e a ~ lb for 146 Gd and the 190 SD bands, and 
(ii) to q — e a+b and p = e a ~ b for the 130 and 150 SD bands; ^ and ^7 are in units of H~\eV. 

Table 2. Theoretical and experimental 7-ray energies and experimental errors for SD bands in 
the A ~ 130 region. Experimental data are taken from Refs. 46-48. 

Table 3. Theoretical and experimental 7-ray energies and experimental errors for SD bands in 
the A ~ 150 region. Experimental data are taken from Refs. 49-53. 

Table 4. Theoretical and experimental 7-ray energies and experimental errors for SD bands in 
the A ~ 190 region. Experimental data are taken from Refs. 44, 45 and 54. 

Table 5. Theoretical and experimental 7-ray energies and experimental errors for SD bands in 
the A ~ 190 region. Experimental data are taken from Refs. 55-58. 

Table 6. Theoretical and experimental 7-ray energies and experimental errors for SD bands in 
the A ~ 190 region. Experimental data are taken from Ref. 59. 

Table 7. Intrinsic electric quadrupole moments for 192 Hg in units of eb. The theoretical moments 
{Qo) q and {Qo) qp are calculated with ^7 = 5.58, (3' = 0.12 10™ 1 and -^ = 5.91, a = —0.15 10~ 2 , 
b = 0.47 10~ 2 for the q- and gp-rotor models, respectively. The experimental values (Qo)ex as 
well as the upper and lower experimental errors AQq and AQq are taken from Ref. 60. 

Table 8. Intrinsic electric quadrupole moments for 194 Pb in units of eb. The theoretical moments 
(Qo) q and (Q ) qp are calculated with ^7 = 5.62, f3' = 0.13 10" 1 and ^ = 5.75, a = -0.78 10~ 3 , 
b = 0.92 10~ 2 for the q- and gp-rotor models, respectively. The experimental values (Qo) cx as 
well as the upper and lower experimental errors AQq and AQq are taken from Ref. 61. 

Table 9. Intrinsic electric quadrupole moments for 196 Pb in units of eb. The theoretical moments 
(Qo)q and (Qo) qp are calculated with 7^7 = 5.68, (3' = 0.11 10 _1 and ^ = 5.71, a = —0.17 10~ 3 , 
b = 0.11 10 _1 for the q- and gp-rotor models, respectively. The experimental values (Qo)ex as 
well as the upper and lower experimental errors AQq and AQq are taken from Ref. 57. 
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Figure captions: 

Fig. 1. The characterization, in the plane of the free parameters a = (3 cos 7 and b = (3 sin 7 of 
the gp-rotor model, of the SD bands in the A ~ 130, 150 and 190 mass regions. 

Fig. 2. Comparison between the theoretical and experimental 7-ray energies in keV. Solid lines 
and dotted lines display the results for the gp-rotor model and the basic model, respectively. 

Fig. 3. Comparison between the theoretical and experimental 7-ray energies in keV. Solid lines 
and dotted lines display the results for the gp-rotor model and the g-rotor model, respectively. 

Fig. 4. Comparison between the theoretical and experimental 7-ray energies in keV. Solid lines 
and dotted lines display the results for the gp-rotor model and the g-rotor model, respectively. 

Fig. 5. Comparison between the theoretical and experimental 7-ray energies in keV. Solid lines 
and dotted lines display the results for the gp-rotor model and the g-rotor model, respectively. 

Fig. 6. Comparison between the theoretical and experimental 7-ray energies in keV. Solid lines 
and dotted lines display the results for the gp-rotor model and the g-rotor model, respectively. 

Fig. 7. Comparison between the theoretical and experimental 7-ray energies in keV. Solid lines 
and dotted lines display the results for the gp-rotor model and the g-rotor model, respectively. 

Fig. 8. The dynamical moments of inertia for 19 °- 192 Hg, 194 Hg(a) and 194 Hg(b) calculated for 
the U q {s\\2) and U qp (u.2) symmetries and compared to the experimental values. The moments 
of inertia are in units of ^MeV" 1 . 

Fig. 9. The geometrical factor G(I) of the reduced transition probabilitiy B(E2) for 192 Hg, 
194 Pb and 196 Pb calculated for the f/ (? (su 2 ), U qp (u 2 ) and su 2 symmetries. 
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